Spin transmission through quantum dots with strong spin-orbit interaction 
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Quantum oscillations of the spin conductance through regular and chaotic 2D quantum dots under 
the varying Rashba spin orbit interaction and at zero magnetic field have been numerically calculated 
by summing up the spin evolution matrices for classical transmitting trajectories. Fourier analysis 
of these oscillations showed power spectra strongly dependent on the dot geometry. For narrow 
rings the spectra are dominated by a single peak in accordance with previous analytic results. In 
other geometries the spectra are represented by multiple peaks for regular QD and quasicontinuum 
for chaotic QD. 
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Recently spin polarized transport has attracted a con- 
siderable interest because of its potential in semicon- 
ductor device applications and quantum computing H. 
In semiconductor quantum wells, such transport can be 
strongly influenced by the spin-orbit interaction (SOI). 
There are two contributions to the SOI, the Dresselhaus 
term Q and the Rashba term j^]. The latter can be 
quite strong in some heterostructures 0]. There have 
been a number of experimental and theoretical studies of 
the spin-orbit effects on the electron transport in semi- 
conductors. Here, we will focus on the quasiclassic elec- 
tron transport through 2D regular and chaotic quantum 
dots (QD). Special attention will be payed to the dou- 
bly connected QD geometries. While considering such 
a problem, an important parameter of the characteristic 
SOI length has to be introduced and compared with the 
QD size L. This length can be defined as Lso ~ 2hvp/A, 
where A is the spin splitting of the electron energy due to 
the SOI and vp is the Fermi velocity. In Grundler's sam- 
ples 0, it has been found as short as lO^^cm. The regime 
of weak SOI can be realized when L <^ Lso- In this 
regime, the random matrix approach has been applied j5| 
to study the SOI effects on the electric transport through 
a chaotic QD 6]. In the opposite case of the strong SOI, 
L > Lso, the spin dynamics in the dot becomes fast and 
the random matrix approach is not helpful. This regime 
is of much importance from the fundamental point of 
view and also for semiconductor device applications, be- 
cause at zero magnetic field it allows to observe periodic 
oscillations of transport parameters due to the interfer- 
ence of spinor amplitudes. Their phases are accumulated 
under the SOI effect during the particle motion along a 
quasiclassic trajectory and depend only on the trajectory 
geometry. In QD having the ring shape constructive and 
destructive interference of the spin phases gives rise to 
periodic variations of the electric conductance as a func- 
tion of the SOI strength 0,0 ■ practice, it is possible 
to change the strength of the Rashba SOI up to 50% of 



its value by tuning the gate voltage Hence, one can 
observe periodic oscillations in the electron transport by 
a simple gate manipulation. 




FIG. 1: (a) 2D electron gas confined to a narrow layer in z- 
direction and bounded inside a two dimensional billiard sam- 
ple in xy-direction is connected to the source and drain by 
two leads. The Rashba effect induces the effective magnetic 
field -ffgg- which is parallel to the xy plane and perpendicular 
to the electron velocity direction. While an electron moves 
along the trajectory, its spin precesses around H^f^- (b) The 
geometry of the samples can be an annular billiard with the 
inner and outer radius ri and r2 and the lead width w or the 
Sinai billiard with radius ri of the inner circle and the length 
2r2 of the outer square. 

For a 2D disordered ring, the spin phase effect on the 
weak localization electric conductance and its pair 
conductance correlation function ^ have been calculated 
under the assumption of a diffusive electron propagation 
inside the ring. Besides the electric transport, the spin 
phase interference also affects the spin transport. Oscilla- 
tions of the mean spin conductance through a classically 
chaotic QD with the shape of a thin 2D loop have been 
analytically calculated [lO|. The width of the loop arms 
was assumed to be much less than Lso- In this case the 
particle oscillations between the loop boundaries produce 
only a small effect on the spin phase. The phase accumu- 
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lates mostly when the particle advances along the loop 
arms making several windings during its escape time. In 
the present work we will calculate the mean spin con- 
ductance beyond the narrow loop geometry. We will em- 
ploy the path integral approach and perform a numerical 
simulation of the spin evolution along the classic trans- 
mission trajectories. Weak localization corrections will 
be ignored under the assumption of a large number of 
transmitting channels. The spin conductance will be cal- 
culated for regular and chaotic systems like the annular 
billiard, the square billiard, and Sinai billiard. 

An appropriate approach to the spin transport through 
a QD sample like that in Fig. ^ is the generalized 
Landauer-Biittiker formula 
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where denotes the transmission amplitude of an elec- 
tron at the Fermi energy propagating from the incoming 
channel n to the outgoing channel m and indices a and 
b stand for the spin variables. The spin conductance 
given by Eq. Q relates the spin current j'^, flowing 
towards the right lead, to the spin polarization matrix 
(j-ac , 7V'A/i/2 in the left lead, where A/i is the difference 
of the chemical potentials between two spin projections 
in the left lead. This difference gives rise to the spin 
polarization along the unit vector TV. Components of 
cr are the Pauli matrices (ct^;, cry, Oz). The transmission 
amplitude can be represented as 

tnni = -i^\/VnVm j dy' j dy(j)*^{y')(t)m{y)Gab{y\y), 

with the transverse wave function (j}.„i {ipn) of the incom- 
ing (outcoming) mode and the retarded Green function 
Gab{y',y), where y {y') denotes a point on the entrance 
(exist) cross section. Under the semiclassical approxima- 
tion, the Green function is represented by a sum over 
classical trajectories propagating from one lead to the 
other Then, according to 10], the contribution of 

such a trajectory 7 to the transmission coefficient can be 
written as 
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where denotes the spin independent amplitude and 
S^'' represents the (a, b) matrix element of the unitary 
spin evolution matrix S-y. In the absence of the external 
magnetic field, this matrix is determined by the Rashba 
spin-orbit interaction and can be written as |lfll | 



(3) 



with the evolution operator S-^. for the i-th straight seg- 
ment of the ballistic trajectory, where S-y-^ starts from 
the entrance lead. For a segment of length k with the 



angle uji with respect to the a;-axis, this operator can be 
represented by the matrix [lCll | 



Sj. = 1 cos ' 
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(4) 

where 1 is the 2x2 identity matrix. The operator Q 
acts on spinors and describes the spin precession around 
the axis determined by the unit vector with coordinates 
(— sinwi, coswi) through the angle 2li/ Lso- The spin ro- 
tation length Lao = h/am* is inversely proportional to 
the Rashba coupling constant a and the effective mass 
TO*. Instead of the 4-th rank spinor, which is the spin 
conductance in the spinor representation in Eq. , it is 
convenient to write it as a 3D matrix 

\ " rahd„ 
9tj ^ Z^^i ^] Sabcd, 

abed 

which determines the spin transport between spin ori- 
entations i and j. In the configurational averaging, the 
spin conductance is averaged over a small range of en- 
ergy around Ep or over small variations of the dot shape 
which smooths out the rapid transmission oscilla- 
tions caused by interference of de Broigle waves. The so 
averaged {gij) is then given by [Toj 



(5) 



One should note a fundamental distinction of the spin 
phase effect, as it shows up in the spin transport, from 
its effect on the charge transport. Indeed, the mean elec- 
tric conductance is given by (O with the trace equal to 
25ij. Hence, it does not depend on the spin-orbit effects. 
On the contrary, the components of the mean spin con- 
ductance oscillate, as we will see below, due to the spin 
phases contained in S^. It is also important to note a 
difference between the AB and spin phase effects. Unlike 
the latter, the AB effect vanishes in both the electric and 
spin mean conductances, as far as the weak localization 
correction is not taken into account. This difference is a 
consequence of non-Abelian nature of the spin phase. 

In the semiclassical approximation, the spin indepen- 
dent transmission probability |t(7)P in ((SJ is given by the 
transmission ratio of an ensemble of classical trajectories, 
in which |i(7)P is replaced by /^cos(0^), where f-y = 1 
or /-y = for transmitted respectively reflected trajecto- 
ries. The value cos(6'7) is the weight of the trajectory 7, 
according to its injection angle 0^, with x axis The 
initial conditions for the trajectories are chosen such that 
Tiw positions are uniformly distributed on the entrance 
cross section with ne uniformly distributed injection an- 
gles between —tt/2 and 7r/2 on each position. The spin 
conductance is calculated as a function of SOI strength 
by setting h = I and taking n^^^ points of Lso between 
Lsoi and Lso2, uniformly distributed on the 1/Lso axis. 
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FIG. 2: The conductance {gyy) as a function of 1/Lao 
and its power spectrum in the ring billiard with 
{r2,w, n^,ng,nL,a , Lsoi , Lso2) = (l, 0.1, 9, 80, 3000, 0.01, 2), 
where ri = 0.9 for (a), (b) and ri = 29/30 for (c), (d). 



Figure |21 (a) shows the conductance {gyy) through 
the narrow ring in the inset of Fig. |2] (b). The con- 
ductance oscillation is regular for small 1/Lso around 
1/Lso < 10 ~ (r2 — ri)^^. In this regime, the spin ro- 
tation length Lso is larger than the ring width and the 
electron spin phase does not change much between two 
collisions with the sample boundary. For 1/Lso > 10, 
the regular oscillation disappears, because the trajecto- 
ries with the same winding number arrive to the exit lead 
with different phases, which generate a disordered inter- 
ference pattern. The power spectrum in Fig. [21 (b) on the 
interval < 1/Lso < 100 shows a dominating frequency 
centered around 0.82 with a low frequency tail. For an 
extremely narrow ring like that in Fig. El{d), the regular 
oscillation regime is extended to around 1/iso = 20. The 
dominating frequency in the power spectrum is shifted to 
0.92. 

For narrow rings the quantum oscillations of the spin 
transmittance can be compared with the analytical result 

m 



(gyy) = 



-.goK^ cos[n ^/T+I{d/lZof] 



(6) 



k2 -^4sin2{cos[7r^l + 4(d/Lso)2]}' 

with the spin independent conductance go, the aver- 
aged radius d — ^ (r\ -I- r|)/2, and the parameter k = 

where r denotes the mean escape time and 
represents the average duration for one winding. The 
equation lO is valid for k and \r2 — ri\/Lso ^ 1- The 
dashed curve in Fig. is rescaled from Fig. |2] (c) and 
is bounded by (gyy) = 0.56, which is the spin indepen- 
dent conductance given by Eq. ^ with tia,ce{aiS^aj SI^} 
substituted for 2. The sohd curve in Fig. 01 is obtained 
from Eq. © with go = 0.56. The positions and magni- 
tudes of the maxima and minima on both curves almost 



coincide for first several oscillations, but the deviation, 
as expected, increases for smaller Lso- Furthermore, the 
dashed curve is much smoother than the solid one, be- 
cause the leads in the ring in Fig. |21 (d) are relatively 
wide and the mean escape time of the particle is not long 
enough to produce the sharp interference pattern from 
many windings. Simulations (not shown in Fig. ^ with 
narrower leads result in curves with the sharp peaks re- 
sembling the analytical solid curve. 
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FIG. 3: The conductance (gyy) calculated from ||SJ with k = 
0.1 (solid) and rescaled from Fig. |21(c) (dashed). 

An apparent trend seen in Fig. |21 from (d) to (b) is 
a faster damping of the dominating periodic oscillations 
with increasing ring width. However, our calculations of 
the spin conductance in annular and circle billiards gave 
an unexpected result, in which new periodic oscillations 
emerge, as shown in the power spectra in Fig. 0] (a) and 
(b) . They are not so prominent as those in the thin rings, 
nevertheless, are rather distinguishable. The power spec- 
trum for the annular billiard in Fig. 0| (a) consists of a 
dominating frequency around 0.32 and several secondary 
peaks. The power spectrum for the circle billiard in Fig- 
ure 0](b) shows that the main peak is still presents near 
0.65. A simple analysis indicates that the trajectories 
contributing to the main peak in Fig. ^ (a) and (b) are 
those with injection angles close to 0. Notably, the oscil- 
lations in the thin rings and those in the annular and the 
circle billiards are fundamentally different. In the former 
case, the oscillations with one dominating frequency are 
of universal nature, because they can be observed in ar- 
bitrary geometry, both in chaotic and regular systems, 
as far as the billiards have the shape of a thin loop ^3 ■ 
However, oscillations of the spin conductance in the cir- 
cle and annular billiards, like in Fig. (a), are entirely 
determined by the specific geometry of these systems. 

Figure 0] (c) shows the power spectrum for the square 
billiard. Unlike the spectra of the ring and annular bil- 
liards, this spectrum consists of a broad range of peaks 
with comparable intensities. Figure 01 (d) shows the 
power spectrum of a Sinai billiard, which is a strongly 
chaotic system. This spectrum is more uniformly dis- 
tributed and less intense in comparison to all other regu- 
lar billiards studied above. It is interesting to note that 
two apparent oscillations of the spin conductance can be 
still observed in this chaotic billiard [Fig. |31(b)], although 
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FIG. 4: The power spectrum in (a) the annular bilhard 
with ri = 0.5 and (b) the circle billiard with ri = 0, for 
(r2, w, n^,ne,nLsa , Lsoi,Lso2) =(1, 0.1, 4, 40, 3000, 0.01, 10); 
and in (c) the square billiard with ri = and (d) Sinai 
bilhard with n — 0.8, for {r2,w,nw,ng,nLs„ , Lsoi, Lso2) 
= {1, 0.2, 4, 100, 20000, 0.001, 100). 



its shape does not look like a thin loop. 




FIG. 5: The conductance oscillation {gyy) of (a) the annular 
billiard from Fig. 4 (a) and (b) the Sinai billiard from Fig. 4 
(d). 



frequencies in the power spectrum. In regular systems, 
the trajectories with close initial conditions stay close to 
each other during their liftime within a billiard and form 
a trajectory bundle, which contributes to strong peaks in 
the power spectra. However, these bundles are absent in 
chaotic systems, because close trajectories exponentially 
fast move away from each other . 
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As seen from above examples, the frequency spectra of 
different billiards reflect their geometry. For the annular 
and circle billiards in Fig. 01(a) and (b), all segments Zj's 
along a trajectory are of equal length [Fig. ^(b)], which 
gives the unique characteristic frequency of all individual 
S'-y.'s in Eq. |0J. When these operators are multiplied in 
Eq. © and |(SJ), a descrete set of the power spectrum har- 
monics is generated. For the square billiard in Fig. 01(c), 
the segments ^i's along a trajectory have different lengths. 
As a result, the broad range of frequencies contribute to 
the power spectrum. In the case of the Sinai billiard, the 
frequency distribution is more uniform compared to the 
spectra of the regular systems in Fig. El(a), (b), and (c). 
Intuitively, this difference can be understood in the fol- 
lowing way. Each trajectory contributes to a finite set of 



